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1. Introduction 

Most of the pioneering enquiries into quantum chaos have been focussed on bounded 
systems — closed billiards, atomic systems — whose classical dynamics knows only a single 
global time scale, the inverse Kolmogorov entropy which describes the ergodic coverage 
of phase space. The absence of other classical times has facilitated the understanding 
of quantum-to-classical relationships and has opened the view for universal features in 
spectra and eigenfunctions of classically chaotic systems. 

Extended systems represent another, still simple, extreme. They reach ergodicity 
only on a time scale that exceeds all other characteristic times. With extended systems, 
basic solid-state concepts enter quantum chaos. The crucial role of long-range spatial 
order, in particular, for spectrum and transport has to be considered in the unfamiliar 
context of dynamical disorder. In a semiclassical approach, some light could recently 
be shed on the spectral signatures of chaotic diffusion, both in the band structure of 
periodic systems |J and in the discrete spectra of disordered systems with localized 
eigenstates 0, |, |. 

A region intermediate between bound and extended is marked by systems comprising 
just a few weakly connected similar cells. In the following, we will consider 
spatial confinements or phase-space structures that decompose into two compartments, 
connected only by a "bottleneck" (Fig. [I]) [[J. Equilibration between the cells then 
takes much longer than the ergodic coverage of a single cell. It constitutes a second 
independent time scale of the classical dynamics. 

Chaotic systems with two cells have quite diverse applications. The phase space of 
alkali atoms in a strong magnetic field may contain two almost disjunct chaotic regions, 
related by a spatial symmetry |J. Similar systems arise as models for non-ergodic 
reaction dynamics in quantum chemistry 0. More generally, they form prototypes of 
systems decomposing into two similar, weakly coupled parts. In this sense, they can 
represent heavy nuclei in a final state of fission ||, or nuclei with their two isospin 
subsystems interacting weakly due to a slight breaking of isospin invariance |J. 




Figure 1. Two-cell billiard with narrow bottleneck, of arbitrary (a) and reflection- 
symmetric shape (b). 

In their quantum-mechanical properties, two-cell systems already exhibit the 
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decisive influence of spatial symmetry. If the two cells are of arbitrary shape or size 
(Fig. [l]a), their restricted communication will merely be reflected in a quantitative 
deviation from the canonical random-matrix statistics ||. It does not introduce any 
qualitatively new feature. The situation changes considerably if the cells are related by 
some twofold spatial symmetry (Fig. |I]b). A genuine quantum phenomenon, a coherent 
mode of transport between the cells emerges, and the spectrum shows systematic 
quasidegeneracies. This clustering of levels is manifest in the two-point correlations 
as a marked positive peak on the scale of the mean single-cell level separation, or in the 
time domain, the single-cell Heisenberg time Jj]]. 

This phenomenon should be carefully distinguished from tunneling. To be sure, the 
doublets do resemble tunnel splittings in the sense that they are based on a discrete 
spatial symmetry and correspond to quantum coherent transport on very long time 
scales. Moreover, in the wavenumber regime where there is no open channel in the 
constriction between the cells, the wavefunctions decay exponentially into this region. 
Even at higher wavenumber, the amplitude is often strongly suppressed there (cf. Fig. |6]). 
Since, however, there is neither a potential nor a dynamical barrier involved, this 
transport is slow but not classically forbidden. 

Chaos-assisted tunneling, in contrast, is a hallmark of bistable systems with a 
mixed phase space [|K], |1^, |TJ|. It occurs between symmetry-related pairs of regular 
islands in phase space that are separated by a chaotic region. Here, by contrast, we are 
dealing with symmetry-related chaotic regions communicating through a narrow bridge 
in space. Still, a similar situation can occur also in mixed bistable systems. In fact, the 
distribution of the splittings of doublet states supported by symmetry-related pairs of 
chaotic regions [T3j| forms an important input to the distribution of tunnel splittings in 



chaos-assisted tunneling |L5 



In the following sections, we shall develop a theory for the spectral statistics of 
two-cell systems that largely rests on recent progress in the analysis of band structures 
of classically chaotic systems with spatial periodicity |TJ. There, an important input 
has been the notion of form factors with a winding- number argument, specific for 
transitions spanning a corresponding number of unit cells. Likewise, the group property 
of a reflection or translation symmetry of a two-cell system enables the definition of 
form factors with a rudimentary spatial resolution, expressed by a binary index that 
indicates either return to the same cell or transport into the opposite cell. 

However, as compared to [Ij], we do not only go from a large value of N, the number 
of unit cells, to iV = 2. Here, we shall concentrate on the case of slow exchange between 
the cells, the "weak-coupling" or "tight-binding" limit in solid-state terminology, while 
IJ was devoted to the opposite case. This implies that, even if the two-cell system is 
"unfolded" to form an infinite chain (see [Appendix A| ), the concept of homogeneous 
diffusion no longer applies. Rather, in the unfolded picture, we are dealing with a 
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spatially discrete diffusion process that deviates significantly, on short time scales, from 
ordinary diffusion. Concerning the spectral statistics, the principal consequence is that 
we are now dealing with flat "bands" . 

In order to extend a semiclassical treatment of the spectral correlations to energy 
scales below the mean level spacing, information on the spreading of the possibly 
complex trajectories that mediate the long-time transitions JL6], [b|, analogous to the 
sum rules for ergodicity |T7) or diffusion [|, would be required. Alternatively, 
we would have to surmount the diagonal approximation. We circumvent this open 
problem and instead adopt a different strategy. In the spirit of Berry's semiclassical 
approximation for the form factors of random-matrix ensembles ||18|| , we impose plausible 
assumptions on the distribution of the narrow splittings. They are based on the relation 
of the doublet splittings in two-cell compounds to the resonance widths in corresponding 
single-cell systems, obtained by opening up the system at the constriction. We derive 
this relation in the case of a single open channel between the cells, where doublet 
splittings and resonance widths are not simply identical. 

Switching back from the distribution of doublet splittings to the corresponding time- 
domain function valid on long time scales, we match the resulting long-time asymptote 
with the semiclassical short-time behaviour. In this way, we achieve expressions for 
the spectral two-point correlations on all time scales. They are universal in that they 
contain, as the only free parameter, the characteristic time for equilibration between 
the cells. It is the two-cell analogue of the conductance, the scaling parameter in the 
case of long chains. 

We shall introduce the classical concepts relevant for the dynamics of systems with 
two identical, connected cells in the subsequent Section In Section [3], we define form 
factors specific for an element of the symmetry group of the system. A semiclassical 
theory for the short-time regime of these form factors is developed in Section |], while 
Section |5] is devoted to their quantum long-time behaviour. Some of these calculations 
are extended to the case of unrestricted values of N in two appendices. Section || serves 
to introduce four illustrative models, two versions of a Sinai billiard [pi], 2~0|j , quantum 
graphs [21] configured in such a way that they form a two-cell system, a two-cell variant 
of the quantum kicked rotor [^, |lj, and a random- matrix model [11, [l2|. Spectral 
data obtained numerically for these models corroborate our theory. Section [7| contains 
a synopsis of the various limiting cases covered in this paper. 



2. Classical dynamics in two-cell systems 

As a minimal version of a classical two-cell system, consider the following model: Two 
spatially confined compartments are connected by some narrow duct (Fig. |I]a). In view 
of the intended applications, we require a few additional properties. The leakage from 
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the cell where the system is prepared (subscript '0' in the following) to the opposite side 
(subscript '1') should be completely described by a single time scale 1/A. This amounts 
to an exponential decay of the population from the initial cell, if it were opened by 
removing the opposite cell. We require the rate A to be the same in both directions. A 
sufficient condition for this to be true is that the cells form a (translation- or reflection- 
) symmetric pair (Fig. [l|b). Finally, we assume that the dynamics within the cells is 
chaotic and thus ergodic, and that coverage of a single cell is reached instantaneously 
on the scale 1/A. 

Under these conditions, the time evolution of the probability to stay in either cell 
obeys the simple pair of master equations, 

P (t) = X(P 1 (t)-P (t)), 

P 1 (t) = A(P (*)-Pi(t)). 

The relaxation into equilibrium of the two probabilities, given by lim^oo P /i(t)/(P (t) + 
Pi(t)) = 1/2, is governed by the rate A = 2A. From an initial state -Po(O) = 1, -Pi(O) = 0, 
they evolve as 

^o/i(t) = ^(l±e- At ). (2) 
The population difference 

Pd(t) = Po(t) - Pi(t) (3) 

is another relevant quantity. Besides the sum Po{t) + Pi{t), it plays the role of an 
eigenmode amplitude of the master equation (1). Its time evolution reads, for the same 
initial state as above, 

Pd(*) = e' At . (4) 

The above considerations apply also if the cells communicate through two or more 
physical channels. This includes in particular the case of two cells connected at both 
"ends" to form a ring. The rates of probability exchange through the channels then 
just add to give the global rate A. The diffusive dynamics that results if a two-cell ring 
configuration is unrolled into an infinite chain, is discussed in Appendix A . 



As an example, we state the explicit expression for the decay rate A in the case of an 
ergodic double billiard as in Fig. |l|. The phase-space area leaving one cell of the billiard 
in time dt through the connecting channel of width s, at unit speed, is dfl = 2sdt. This 
is to be normalized by the area A of a cell and by 2tt, the size of momentum space 
projected onto the energy shell. The resulting expression for the escape rate is 
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3. Generalized form factors 

In quantum systems, each unitary symmetry gives rise to a constant of the motion, 
a "good quantum number". It takes as many values as there are irreducible represen- 
tations of the symmetry, and the full spectrum can be decomposed into subspectra, each 
of which pertains to a given irreducible representation. Formally, the decomposition 
is effected by the projectors |23| P u = N^ 1 J2n=o Xv{9n)U\g n ), v — 0, . . ., N — 1. 



Here, N is the number of elements g n of the symmetry group and simultaneously the 
number of its representations (for simplicity, we assume all representations to be one 
dimensional). The character of g n in the z/th representation is referred to as Xv{gn), and 
U(g n ) denotes the unitary transformation corresponding to g n . Spectral densities and 
correlation functions within a given representation can then be defined on basis of the 
symmetry-projected Green function G U (E) = P U G(E), with G(E), the Green function 
for the entire spectrum. For example, the symmetry-projected spectral density in the 
vth representation is defined as 

1(E) =J2 5 ( E ~ E *,») = --Imtr [G V (E)\, (6) 

71 

where the E a ^ u are the eigenenergies in the ^th representation. 

For the study of quasidegenerate doublets, another type of spectral density is 
as relevant as the symmetry-projected one. As an alternative to d u (E), one may 
define densities and derived quantities that refer to a group element g n , instead of an 
irreducible representation v. The symmetry group induces a tiling of (phase) space, i.e., 
a decomposition into disjunct segments such that each of them is mapped onto all the 
others by the transformations in the group, thus covering the entire space ||24j| . A group- 
element-specific spectral density therefore provides a rudimentary spatial resolution on 
the scale of the fundamental domain of the group. 

If all representations are one dimensional, the set of columns of the matrix \v,n = 



Xu(gn) of group characters forms an orthogonal basis in N dimensions [23]. The same 
is true for the rows. The matrix as a whole therefore has full rank and is invertible. We 
refer to the inverse matrix as x~ l - Multiplying the vector of symmetry-reduced spectral 
densities d u (E), Eq. (|J), from the left by x~ X i we obtain the spectral densities ]T] 

N-l 

d n (E)= Y.(x- 1 Undu{E) (7) 
= --Im / dqG(g n (q),q;E). (8) 

71 Jfd 

The space integral in Eq. (8) extends only over the fundamental domain (subscript 'fd') 
of the tesselated space. The first argument of the Green function in the third line is the 
image of its second argument, q, under g n . This suggests that d n (E) refers to transitions 
from any segment to its image under g n . 
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By Fourier transforming Eq. (7) with respect to energy, we arrive at the analogous 
density in the time domain, 

/oo 
dre-^-^r/Hd)) (9) 
-oo 

= / dq(g n (q)\U{tnT)\q). (10) 

We have switched to dimensionless energy, r = (d {d )E, and time, r = t/tu, by 
scaling with the mean spectral density (<ifd) in the symmetry-reduced space, and the 
corresponding Heisenberg time tu = 2iih(dfd), respectively. Equation (10) describes an 
amplitude to return modulo the symmetry transformation g n . The corresponding return 
probability is given by the form factor 

K n (r) = ^-\a n (r)\ 2 , (11) 

where Ar^ is the total width, in units of (rffd), of the spectrum considered. Eq. ([□]) 
is equivalent to a definition of the form factor as the Fourier transform of a spectral 
autocorrelation function |Q [|. 

For a two-cell system, a symmetry that maps one cell to the other can be a reflection 
or a translation. Their group elements are identity (denoted by n = in the following) 
and reflection or translation (n = 1, without distinguishing the two). The characters are 
both 1 in the symmetric (subscript '+') representation, and ±1 in the antisymmetric 
('— ') representation. Following the general discussion above, we define symmetric and 
antisymmetric form factors, respectively, by 

K ± {t) = -^-\d ± (r)\\ (12) 

The amplitudes a±(r) are obtained, for example, by sorting the spectral data of the 
two-cell system according to the symmetry of the corresponding eigenstates, and Fourier 
transforming as in Eq. (9). Alternatively, K± can be interpreted as the form factors 
of a single cell with Neumann or Dirichlet boundary conditions, respectively, imposed 
on the line(s) in configuration space common to the two cells. In general, this choice 
of boundary conditions does not affect the chaoticity of the classical dynamics. Hence 
we expect the form factor of a single cell to equal the random-matrix result ]25] up to 
normalization, i.e., 

K±(t) = ^rmt(t) . (13) 

Form factors specific for return to the initial (subscript '0') or switching to the opposite 
cell ('1') are defined according to Eq. (7) and Eqs. (9), ( |TT| ) as 

K 0/1 (T) = ^-\~a + (T)±~a_(r)\ 2 . (14) 
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Note that the amplitudes are superposed before squaring. We shall show in Sect. ^| 
below that considering the incoherent superpositions K (t) ± K 1 (t), in turn, provides 
an approximate access to the distribution of doublet splittings and inter-doublet 
separations, respectively. 

The combination of amplitudes a + (r) + a^(r) entering Kq(t) in Eq. ( [H]) is obtained 
in (9), when d n on the r. h. s. is replaced by the total spectral density of the two-cell 
system. Hence K Q is — up to scaling of time and energy — equivalent to a form factor 
defined without any reference to spatial symmetry of the two cells. 

In Eqs. fljjD and ( |l4t ) we have chosen a normalization which ensures that K±{r) 
and Kq/i(t) approach the same value 1/2 for r — > oo, provided the two cells are not 
completely disconnected. In addition, the sums of the symmetry-projected and the 
group-element- specific form factors are the same for arbitrary time r, 

K (t) + K x {t) = K+(r) + K-{t) . (15) 

This identity may be interpreted as a preservation of norm and follows generally from 



the unitarity of the matrix \v,n of group characters. Because of fll3|) it leads to the 
relation 

K (r) + K^r) = K RMT (r) . (16) 

It is instructive at this stage, to consider the trivial limiting cases of two completely 
isolated or two very strongly interacting cells, respectively. In the first case, we have 
from the definitions (7), (9) and flU) K x (t) = 0. Then (0) implies K {t) = K rmt {t), 
and this is indeed what is expected within our scaling of time and energy from the fact 
that the total spectrum is the superposition of two identical random-matrix spectra. In 
the other extreme, the two subspectra of positive and negative parity can be considered 
statistically independent, which has the consequence K (r) = K\{t) = K rmt (t)/2. 

4. Semiclassical regime 

In order to construct a semiclassical trace formula for the symmetry-projected spectral 



density d„(E), Eq. (y), the concept of periodic orbits has to be extended |2J]. In case 
the dynamics within the cells has no significant admixture of regular motion, the trace 
formula reads 

1 t-i(p) 

x X u{gj) exp - ifij^j (17) 

The sum now runs over generalized period orbits j. Their end point is not necessarily 
identical with the starting point, but must be mapped to it by some group element 
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Qj. The corresponding term in Eq. (17) then contains the character Xv{9j) as an extra, 
non-classical phase factor. A correction of the amplitude for orbits that coincide with 
symmetry lines is effected by Kj As usual, Tj , Mj, Sj, fij denote primitive period, 
stability matrix, classical action, and Maslov index, respectively, of orbit j. 

The role of the generalized periodic orbits becomes even more transparent in the 
analogous trace formula for the group-element-specific density, 

1 rp(p) 



mN i <W| det(Mi - /)| 



/ S 7T \ 

x %j, 9n) exp M-^ - ifXj-j . (18) 

The delta function in the second line equals unity if its arguments coincide and vanishes 
otherwise. It selects orbits j whose endpoints are connected by g n . They mediate 
transport from the original segment to its image, with the restriction that initial and 
final points are exactly related by the symmetry. 

The interpretation that spectral quantities associated with g n describe transport 
from an original space segment to its image under g n is borne out quite explicitly by the 
form factors. A semiclassical expression for the K n (r) can be derived by substituting into 
Eq. flTT|) the trace formula (18), Fourier transformed to the time domain as in Eq. (9). 
Within the diagonal approximation with respect to pairs of generalized periodic orbits 
|Tj| [I], fjj, which is valid for times t -C tu_, one obtains, 

K { : c \t) = ln TP(g n , Tt u ), r « 1. (19) 

Equation ( |I9"D relates the form factors to the classical probability P(g n ,t) to return 
in time t to a phase-space point related to the starting point by g n . The contribution 
of repetitions of shorter periodic orbits has been neglected in Eq. ([19]) • By introducing 
a global degeneracy factor 7„ to account for antiunitary symmetries like time-reversal 
invariance, we ignored the occurrence of self-retracing orbits. This factor takes the value 
2 if orbits that are periodic modulo g n are generically time-reversal degenerate, and 1 
else. A non-trivial dependence of 7 n on n can occur, e.g., in periodic systems with 
iV > 3 unit cells |l|. There, time-reversal invariance is generally broken for orbits with 
winding numbers nmodiV ^ 0, N/2, due to Bloch phases that are not real. 

In the spirit of the known classical sum rules for ergodic systems ||17|| , we assume that 



the generalized periodic orbits are not distinct from the generic non-periodic ones in their 
average spreading. We can then relate the P(g n , t) to the classical propagator p(r', r; t) 
(the integral kernel of the Frobenius- Perron operator) by a phase-space integration over 
the fundamental domain, 

P{g n ,t) = I drp(g n (r),r;t), (20) 

Jfd 
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where r = (p, q) denotes a phase-space point within the fundamental domain on the 
energy shell. In case that the chaotic coverage of the single cells is homogeneous, the 
classical propagator depends on g n but not on r and can thus be expressed by the coarse- 
grained probabilities defined in section 0. In this case the integration in Eq. ( p0|) 
becomes trivial and results in 

P(g n ,t)=P n (t). (21) 



For the group-element-specific form factors defined for two-cell systems, Eq. ([14]) , the 
semiclassical expression finally reads 

K^(r)= 1 rP 0/1 (rt u ) (22) 

= ^(l±e- 2A *^), r<l. (23) 

We have used here that in symmetric two-cell systems, the degeneracy factor 7 does not 
depend on g n . 

A detailed investigation of the phase-space coverage for a specific double billiard 
p8fl shows that — for a finite time depending on the special properties of the employed 
model — there can be small deviations from the homogeneity assumed in the derivation 
of Eq. (22). Similar restrictions of our theory may arise from the neglect of marginally 
stable (bouncing-ball) orbits. However, all the approximations discussed so far are 
standard within a semiclassical theory for two-point correlations, and although they 
cannot be rigorously justified, they are sufficient to reproduce most of the available 
numerical data [0, [|, |, § §. 

The most important limitation in this respect is due to the diagonal approximation 
for systems with time-reversal invariance, 7 = 2. It is well known that this approxi- 
mation reproduces only the slope near r = for systems with a single chaotic cell, and 
it is not surprising that we observe the same for a two-cell system when we compare (23) 
to Eq. ([16]). The semiclassical result deviates exactly by the same factor 7r/i^ RM T( r ) 
known from simple ergodic systems . 

In the remaining paragraphs of this section, we will discuss briefly the influence of a 
small breaking of the symmetry of the two billiard cells. It is clear that the overshoot of 
K (t) over its asymptotic value, as described by Eq. (23), is the result of the symmetry. 
For an asymmetric double billiard, this overshoot should vanish altogether, leaving the 
steeper initial rise of the form factor implied by Eq. (23) as the only spectral signature 
of the restricted exchange between the cells. However, we expect that there exists a 
continuous crossover, as a function of some parameter, that expresses the degree of 
symmetry breaking. 

A semiclassical approach that extends the above arguments to the case of a 
weakly broken symmetry has been presented in [|S5[| . It is based on the idea that the 
contributions of a given set of N symmetry-related periodic orbits to the form factor will 
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no longer be N-fold degenerate, but can still be completely included, without resorting 
to a diagonal approximation within this set. The following assumptions have to be 
made to justify this strategy: (i), the perturbation is sufficiently weak not to destroy 
the structural stability of the periodic orbits, i.e., no periodic orbits disappear or are 
created, as compared to the unperturbed system; (ii), only the change AS in action 
has to be taken into account because it appears in the exponential, while the changes 
in amplitude and period can be neglected; (iii), the changes in action result from many 
statistically independent perturbations of the orbit so that, by the central limit theorem, 
they can be considered as Gaussian random variables with zero mean and variance 

{(AS) 2 } = 5 2 r. (24) 

The constant 5 will serve as the basic parameter for the degree of symmetry breaking. 
The proportionality to time reflects the accumulation of squared action changes over 
the length of the periodic orbit. 

Under these assumptions, the influence of symmetry breaking can be expressed as an 
effective degeneracy relating the form factor for broken symmetry to the corresponding 
one for the symmetric system (in the present context, we are left with the "symmetry- 
insensitive" Kq(t)). For a two-cell system, we find 

Kt ym) (r,5) = ^l + e-^)K (r). (25) 
Inserting Kq C \t) from Eq. (23), this implies 

Kt 7m \r, S) = hl + e- 52r )(l + e- 2A ' HT )- (26) 

Since this result involves the same semiclassical approximations as we made in the 
derivation of Kq(t), its validity is likewise restricted to the short-time regime r < 1. 
Within this regime, it provides the seaked interpolation between the symmetric limit, 
5 2 <C 1, and the limit of totally broken symmetry, <5 2 <;1. Due to the assumptions 
enumerated above, in particular that of structural stability, we expect Eq. ( p5| ) to become 
unreliable also in the latter limit. 



5. Long-time regime 

In the semiclassical time range, we succeeded to express the form factors for all 
parameter regimes, from two-cell systems without significant separation of the cells 
to pairs of nearly uncoupled cells, by a single expression, Eq. (23). We cannot achieve 
this generality for the regime of long times t ^ £h- The case of weakly coupled (Xtu -C 1) 
symmetric billiards, to be considered here, requires additional input besides the classical 
information contained in Eq. (23). At the same time, this is the most interesting 
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situation because only here, quasidegenerate doublets occur with a splitting much 
smaller than their typical separation. 

As a starting point for an alternative approach valid in the long-time regime, 
we return to the exact definition of the group-element-specific amplitudes ao/i(r). 
Specializing Eq. (9) to the two-cell case and inserting the definition (^) of d / 1 (E), 
we obtain 

i i N d 

«n(r) = -E e7riW E e " 27rirrQ ^ n = 0,l. (27) 

^ u =0 a=l 

Here, is the number of doublets in the spectrum, i.e., half the total number of levels. 
We have used the fact that the inverse characters for the twofold reflection or translation 
group can be concisely written as (x _1 )^n = e nmu , with v = 0, 1 corresponding to the 
symmetric and the antisymmetric representations, respectively. In these representations, 
returning to the symbols '+' and '— ', r a $/i = r a; ± = (d id )E ai ± denote the scaled 
eigenenergies. 

We introduce the concept of doublets by writing the eigenenergies as 

r Q ,± = R a ±r a . (28) 

The long-time limit of the form factors for the canonical random-matrix ensembles 
is usually derived under the assumption that the full phases rr are random for 
t ^> 1. Likewise, we here assume the analogous phases rR a contributed by the doublet 
midpoints to be random in the long-time limit. Upon squaring the amplitudes a±(r) to 
obtain the corresponding form factors, this amounts to a diagonal approximation with 
respect to the index a, 

1 1 Nd 

K /i(r) = ^±xr r E cos(47rrr a ). (29) 
In fact, Eq. ( p9|) can be derived also if the two-cell system is unrolled to an infinite chain 



and the doublets are considered as points of continuous bands, see |Appendix B . 

If iVd is sufficiently large, we can replace the sum over a by an integral and consider 
the integration as a Fourier transformation, to obtain 



ffo/i(T) = i(l±j>d(2r)), (30) 



with 



oc 



Pd( T ) — I dr cos(27rrr)pd(^)- (31) 



o 



This is the Fourier transform of the distribution of doublet splittings 

Pd(r) = E <K r - \ r a - - r a<+ \) . (32) 

a=l 
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Forming the difference of the two form factors, 

Pd (r) = K (r/2)-K 1 (r/2), (33) 

we see that this equals the time-domain splitting distribution in its long-time, or 
equivalently, low-energy limit. Given that the form factors contain information merely 
on two-point correlations irrespective of symmetry, it is actually surprising that they 
can be related, as in Eq. ([33|), to a quantity that requires an unambiguous identification 
of doublets. This can be explained by the fact that we had to assume in the derivation 
that the midpoints R a of the doublets are statistically independent of their splittings r a , 
which requires a clear separation of scales between splittings and spacings of doublets. 
Indeed, Eq. ([33]) ceases to be valid for r ^ 1, corresponding to the regime of large r a ~ 1. 

We do not have any semi classical access to Pd{r). In order to nevertheless make 
some progress, we shall resort to results of random-matrix theory on the distribution of 
resonance widths, and argue that the doublet splittings obey a similar distribution. 




Figure 2. Two identical scattering systems are obtained when the channel connecting 
the two halves of a symmetric two-cell billiard is replaced by a semi-infinite waveguide 
of constant width. We argue that small doublet splittings between states of the two- 
cell billiard correspond to narrow resonances of the corresponding pair of scattering 
systems (see text). 

Suppose the channel between the two cells to be replaced by a semi-infinite duct of 
constant width, so that two single cells remain, each with a small opening that couples 
its interior to the continuum in the duct. This situation is illustrated with a symmetric 
two-cell billiard in Fig. |2|. The spectra of the open single cells will then exhibit narrow 
resonances at roughly the same energies where the corresponding closed two-cell system 
shows doublets. It is plausible that the doublet splittings of the two-cell configuration 
are related to the resonance widths in the single-cell setup, at least in a statistical 
sense. Indeed, this assumption is often made, e.g., in nuclear theory, and supported by 
semiclassical and random-matrix arguments. In short, it is justified by the fact that both 
quantities, widths and splittings, can be expressed by the same wavefunction overlaps 
and therefore should obey the same distribution. 
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Here we are interested in the case of very few open channels in the connecting 
section, since this is required for the formation of narrow doublets. In this limit, we 
cannot expect random-matrix theory to remain valid. However, taking the point of view 



of the scattering approach to quantization |Pq| , we show in the following paragraphs 
that there is still a close relation, though not an exact identity, between the respective 
distributions of widths and splittings. 

We shall consider specifically the wavenumber regime where there is just a single 
open channel in the connecting section. This case is realized in the majority of the 
numerical examples below. Here, the calculation is particularly straightforward and 
transparent, because the scattering in the cell is described by a 1 x 1 scattering "matrix" 
(since there is only one opening, all incoming waves leave as reflected waves). As we 
are interested in isolated narrow resonances, we assume their widths to be small as 
compared to their separations. The S matrix can then be written in the form 

S(r) = ^ r - r °- i9 , (34) 

where ro — ig denotes the position of the resonance pole in the complex energy plane, 
again in units of the mean level spacing. The quantization condition in terms of S(r) is 



3(r ± ) = ±1, (35) 

the upper sign referring to the case of symmetric states (subscript '+'), the lower one 
to antisymmetric states (' — ')• Inserting Eq. ( |3"4"D gives the eigenvalues 

r + (9) =r + g cot 9/2, (36) 

r_(0) = r -gtan9/2, (37) 

for the symmetric and antisymmetric cases, respectively, separated by the splitting 

r(0) = \r-(0) - r+{6)\ = . (38) 

| sin 9 | 

This function is it periodic. In the interval < 9 < n, it has a minimum at 9 = 7r/2, 
with functional value r(jr/2) = 2g. It diverges at 9 = 0, tt. 

This analysis already exhibits the essential facts to be demonstrated: There is a 
connection between resonance widths and doublet splittings, but it depends on the 
unknown value of the total phase 9 of the S matrix. As no value of 9 is singled out a 
priori, we assume equidistribution of the total phase. Under this condition, the main 
contribution to the distribution of the splittings comes from r s ^ 2g. This is the simple 
relation between doublet splittings and resonance widths we seek. Quantitatively, we 
find the probability density 

Pd(r|s) = ^(l-(^) 2 ) ' (r>2g). (39) 
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for the splittings. It is normalized to unity, but already its first moment diverges. 
Indeed, as we started from the assumption of small splittings, we cannot expect the 
result to be valid for large splittings. The missing cutoff will be given by our semiclassical 
considerations which cover the complementary regime of large splittings. 

The distributions of wave-function amplitudes and of resonance widths are among 
the established principal results of random- matrix theory |27j. Even if details of their 



application to quantum chaotic scattering are still under study, we can, for the present 
purposes, adopt the canonical random- matrix results for p(g), the probability density 
of the resonance widths, and substitute them to obtain the unconditional splitting 
distribution Pd( r )- 

Its general relation with the conditional Pd{r\g) and p(g) reads 

POO 

Pd(r) = / dgp(g)p d (r\g). (40) 
Jo 

In order to get back from the energy to the time domain, we perform a Fourier 
transformation of Pd(r), cf. Eq. (|31~1) . Inserting the explicit expression (|3~^) for Pd{r\g), 
we obtain 

If 1 2 /2r\ 
p d (r) = - dx r p — , (41) 

7T JO VI - X 2 \ X J 

where p(r) = J °° dg cos(2ir gr)p(g), in turn, is the Fourier transform of the distribution 
of resonance widths. 

For time-reversal-invariant systems, the resonance widths obey a Porter-Thomas 



distribution |27 



e -9/2<9) 



p PT (g) = -j=— ■ (42) 
'2n(g)g 



(43) 



The integral obtained by inserting the Fourier transform 



Ppt(t) 



l + ^l + [An(g) T 

\l 2(1 + MsH 2 ) 



2 



into (|4l|) can be used for a numerical computation of the form factor. In the long-time 
limit we find from an asymptotic expansion of (43) 

If time-reversal invariance is broken, the resonance widths are exponentially distributed, 

Pe xp G?) = ^ye-^\ (45) 
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Pcxp(r) 

and accordingly 

Pd(r) = 



1 + (271(g), 



1 + 47r(<7>T^l + (47r(<7>T)2 + (A-K{g)rf 



(46) 



(47) 



For (g)r > 1 we find p d (r) -> (47r(#)r)~ 2 /2. 

Note that for both, exponential and Porter-Thomas distribution, the asymptotic 
behaviour of the doublet splittings for large time/small energy is equivalent to the 
corresponding resonance- width distribution up to a constant prefactor which relates the 
mean doublet splitting to (g). This constant — it equals 27rr 2 (|)/r 2 (|) « 3.44 with and 
Vo ~ 2.82 without time-reversal invariance — is somewhat above 2 as anticipated from 
Eq. (H). 



0.1 




0.1 
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Figure 3. Time-domain splitting distribution Eq. ( |4l|) for symmetric two-cell systems 
with a single open channel in the connector, in the presence of time-reversal invariance 
(upper solid line) and in its absence (lower solid line). For comparison, the broken 
lines show the Fourier-transforms of the corresponding resonance distributions, i.e., of 
the Porter-Thomas distribution Eq. (43), above, and of the exponential distribution 
Eq. ([46|), below. The parameter value common to all curves is (g) =0.1. 

In Fig. |3], we compare the Fourier-transformed unconditional splitting distribution 
Eq. ([11]) in the presence and absence of time-reversal invariance to the corresponding 
distribution of resonance widths for (g) = 0.1. We see that for (g)r^l, the deviation 
between the two is not dramatic and conclude that the resonance distribution p(g), 
Porter-Thomas or exponential, is the crucial input for Pd(i~), while it is quite robust 
against changes and approximations entering via p&(r\g). 

Now we return to our main line of reasoning and attempt a matching of the short- 
time (large-separation) with the long-time (small-splitting) regime of the form factors. 
This will simultaneously allow us to calibrate the as yet undetermined parameter (g) of 
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the resonance distribution with respect to the classical decay rate A. We shall present 
this calculation only for the simpler case of broken time-reversal invariance. If time- 
reversal symmetry is obeyed, the bad performance of the diagonal approximation at the 
Heisenberg time makes an analogous procedure more problematic. We will discuss this 
in connection with our numerical results in Sections p.'2\ and |6.3j . 

From the semiclassical side, Eq. (23), we find at the matching point r = 1, 

^o/i(l) = ^(l±e- 2A ' H ), (48) 



while from the long-time side, substituting Eq. ( i?D in Eq. (p0[), we have 



K 0/1 (l) = Ul± - ; y=L= - ) , (49) 



2\ l + cVTT^ + c 2 , 

introducing the shorthand c = 8ir(g). These equations are consistent with one another 
if 

e -2Xt H = 1 (5Q) 

or, resolving for c, 

1 _ e -2At H 



2Xt H (0 _ p-2Ai H ' 



^/e- 2A *H(2 



(51) 




(52) 



For A£h 1 (narrow connecting channel or weak coupling), proportionality (g) ~ 
A£h/47t results. It represents a simple relation between the parameter of the quantum- 
mechanical splitting distribution and the classical time scale of equilibration between 
the cells. 

We state the full time dependence of the form factors in the long-time regime, again 
using the abbreviation c for the sake of conciseness, 

i ) 

1 + crVl + c 2 t 2 + c 2 r 2 J ' 

In Fig. £|, we give a synopsis of -K" /i( r ) ( a ) an d ^o/i( r ) = /o°° ^ r [1 ~ -K"o/i( r )] cos(27rrr) 
(b), for values of the decay constant ranging from A <C 1 to A « 1. The figure illustrates 
the crossover of the spectral two-point correlations from the regime of almost immediate 
equidistribution between the cells, A^l, where the two-point statistics barely deviates 
from the corresponding GOE or GUE prediction (for the figure, we have chosen the 
case of broken time-reversal invariance where the semiclassical approximation to the 
random- matrix form factor is exact), to the regime of weak coupling, A 1, with Kq{t) 
rising to a marked peak near r = 1. 
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Figure 4. Group-element-specific form factors (a), as described by Eqs. (23) and (|52|), 
and corresponding cluster functions (b), for the case of broken time-reversal invariance. 
In (a), the upper curves show Kq(t), the lower curves K\{t). In (b), the graphs with 
positive initial slope correspond to Yq(t), those with negative initial slope to Y\(r). 
From the outmost to the innermost pair of curves, the decay rate takes the values 
Atn = 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1.0. Graphs of form factors and cluster functions 
for equal values of Xtu share a common line signature. 



6. Models and numerical results 

In the following subsections, we shall introduce five quite diverse models that allow to 
construct systems with two coupled compartments. The numerical results obtained for 
these models serve to illustrate and check various aspects of the theory developed above. 

6.1. The Z-shaped billiard 




H 
1 



Figure 5. The Z-shaped billiard and the three different families of bouncing-ball 
orbits. 

We construct a two-cell billiard from two quarters of a Sinai billiard |19[ and a 
straight channel such that the resulting shape resembles the letter Z [^8] (Fig. |^). The 
width of this channel will serve as the basic length unit. The remaining parameters 
of the billiard are then the length L of the channel and the common radius R of the 
quarter-circle sections of the boundary. 
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Figure 6. Contour plot of the absolute square of a pair of eigenfunctions of the double 
billiard with quasidegenerate energies. The state with even symmetry has k = 3.7569 
(top), the antisymmetric state has k = 3.7576 (bottom). The geometric parameter 
values are R = 10 and L = 5. 



Since the billiard boundary consists exclusively of defocussing and neutral 



components, the classical dynamics is ergodic and mixing [15]. Hence we can assume 
that Eq. (0) holds to a good approximation, although for finite time systematic 
deviations from ergodicity, e.g., due to the presence of bouncing-ball orbits (Fig. ffy, 
can be observed. In the following, we neglect such effects which have been studied in 
detail in ||. 

The employed quantization scheme is described in |Appendix Q . Figure |] shows 
a representative example of a pair of eigenfunctions of the double billiard with 
quasidegenerate eigenenergies. 

Figure [7| shows a comparison of the various form factors, defined in Section 
with random-matrix theory. In Fig. |7|, we show the symmetry-projected form factors 
K + {r) (panels a,b) and K_{t) (c,d) (cf. Eq. ([12])), as well as the group-element-specific 
ones K (t) (panels e,f) and Ki{t) (g,h) (cf. Eq. (|14D), together with the corresponding 
predictions of random-matrix theory. 

The right-hand column contains blowups of the short-time regime of the data shown 
on the left. The K±(t) follow closely the GOE prediction for a single cell of the double 
billiard, as expected. For short times, Kq(t) and Ki(t) deviate significantly from the 
GOE shape, in the way predicted by semiclassical considerations. The initial slope of 
Kq(t) is increased by a factor 2, while that of K\{t) vanishes. However, the amount of 
data obtained for this model is too small to allow for a quantitative comparison with the 
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Figure 7. Time evolution of various form factors (heavy lines), compared to random- 
matrix theory (broken). The upper two rows show the symmetry-projected form factors 
(panels a - d) , the lower ones the group-element-specific form factors (e - h) . The right 
column consists of blowups of the short-time regime of the data shown on the left. The 
geometric parameter values are R = 5 and (L) = 6.7. 

semiclassical theory beyond the vicinity of r = 0. Moreover our numerical quantization 
procedure did not allow to go to parameter values where Atn~0.5 such that K (r) 
is expected to overshoot near r = 1. For the parameter values underlying the data 
shown, we have Xtu « 4, so that classical equilibration between the cells occurs around 
r w 0.125, cf. Eq. (23). 

6.2. The Sinai billiard 

Another example of a two-cell billiard is provided by one half of the Sinai billiard as 
shown in Fig. |8|. A semicircle of radius R divides a rectangle with side lengths L x , 
L y , into two parts connected by an opening of size s = L y — R along the symmetry 
axis. According to Eq. (|5|), the width of this constriction determines the classical rate 
of transitions between the two cells in the ergodic regime. The geometry of this billiard 
differs from that of the system discussed in the previous subsection in three respects: The 
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single cells are no longer symmetric in themselves, there is no extra connecting channel of 
variable length, and the full configuration has reflection rather than inversion symmetry. 
The main advantage is, however, that there exists a more efficient quantization algorithm 
f20|j , again based on the scattering approach ||26|| . The reflection symmetry allows to 
compute the eigenvalues in the two parity classes separately by requiring Neumann or 
Dirichlet boundary conditions along the symmetry axis. We unfold both spectra using 
the area and circumference contributions to the mean spectral density of one cell 
and arrive at the scaled energy eigenvalues. 

Our theory, developed in Sections @ to H is based on spectral two-point correlations 
that indiscriminately include all level pairs in the spectrum. The symmetry-based 
quantization procedure used for the present model, by contrast, gives us immediate 
access to the scaled eigenvalues r a> ±, presorted according to parity. We take this 
opportunity to make a few remarks concerning the "genuine" doublets, i.e., level 
pairs with identical quantum number a but opposite symmetry, and their splittings 
r a = r a _ — r a)+ . We emphasize again that only in the regime of small splittings, 
statistically independent of the positions of the doublet centers, the two-point statistics 
embodied in the form factors coincides with the distribution of the genuine doublet 
splittings, cf. Eq. (p2|). Outside this regime, the two-point statistics includes separations 
that are possibly very small but belong to states labeled by different quantum numbers, 
and therefore do not contribute to the splitting distribution. In effect, the two-point 
statistics is less restrictive and shows more weight at small separations than the splitting 
distribution. 




Figure 8. One half of the Sinai billiard consisting of a rectangle and an inscribed 
semicircle which divides the system into a reflection-symmetric pair of cells. The 
numerical data presented in this section correspond to L x = 2, L y = 1, and s — 0.05. 

Figure || shows the individual doublet splittings (dots) and a running average (wiggly 
line) as a function of the energy. We observe that the average splitting essentially 
depends on the number A = [fes/V] of open quantum channels in the constriction. For 
low energy, quasidegenerate doublets prevail. In particular, below the threshold energy 
of the first quantum channel, we have \r a \ < 1 for all pairs of eigenvalues. Because of the 
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Figure 9. Spacings between pairs of unfolded energy eigenvalues with equal quantum 
number with respect to a single cell in Fig. || and Dirichlet/Neumann boundary 
conditions along the symmetry axis. The vertical lines mark the threshold energies 
where the first three quantum channels between the cells open. A running average 
of the spacings (wiggly line) is compared to the high-energy approximation (smooth 
curve). 



analogy with actual tunnel splittings |L4j], we presume that a semiclassical description 
of the spectrum in this regime should include also orbits with complex action including 



the diffractive orbits studied in pOfl . This question will be investigated elsewhere. 

As the energy approaches the opening of the first channel, the mean doublet 
splitting increases exponentially, and doublet splittings larger than the mean level 
spacing accumulate. Beyond the opening of the second quantum channel, even the 
average splitting exceeds the mean spacing. Consequently, for high energy, the notion 
of doublets becomes irrelevant for the spectral statistics of the composite system. It is, 
though, well suited in the regime of, e.g., a single open quantum channel, as we will 
show below. 

An approximation to the mean value of the doublet spacing is obtained from the 
asymptotic expansion of the mean spectral staircases N± of the two subspectra. While 
the leading contribution depending on the area A is the same for both spectra, the second 
term depends on the circumference u and the boundary conditions. For h = 2m = 1, 
we have 

A ' ' U -E 1 ' 2 NJE) = NJE) - —E 1 / 2 . 



N-(E) 



-E 



An 4tt 2tt 
With the approximate quantization condition for scaled energy, N±(E ai ±) 
3l||, this leads to 



(53) 
a + 1/2 

(54) 



which is represented by a smooth solid line in Fig. [5| We see that for low energy, the 
approximation (Ell) is correct only in the vicinity of the channel openings, while the 
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mean splitting is approximately constant between the thresholds. Accordingly, taking 
the value predicted by Eq. (|54|) at the opening of the first channel r = An /4s 2 for the 
entire subsequent interval till the next threshold, we find that the mean dimensionless 
doublet splitting for one open channel is 1/2, independently of the size of the hole. Thus 
it is already of the order of the mean level spacing of the composite two-cell system. 

It is an important point that this does not restrict the applicability of our theory: 
The high probability of large doublet splittings corresponds to the fact that Pd( r )> as 
obtained from Eq. (|3"PD, has a diverging first moment. Nevertheless, its Fourier transform 
is well behaved. Beyond the Heisenberg time, where we make use of it, it is essentially 
determined by the behaviour of the distribution at small spacings. Indeed, Fig. |9] shows a 
large number of doublets with a width well below the mean level spacing, which justifies 
our approach. 




Figure 10. (a) Form factors i^ /i( T ) f° r the two-cell Sinai billiard of Fig. |[ The 
smooth solid curves represent the semiclassical result in diagonal approximation (23), 
the dashed lines show a fit to the ansatz which explicitly obeys Eq. (|T^). The 
validity of the sum rule ( |l6|) is demonstrated in (b), where the solid line represents the 
GOE form factor. 



In Fig. [iO|a, we present the form factors K /i(r) obtained from the 1,187 doublets 
with one open quantum channel. They were computed after splitting the spectrum into 
small intervals of 30 doublets each. For the parameter Atn entering the semiclassical 
theory we use the value obtained from Eq. (|5|) with k as at the center of the 
considered interval. The dashed line shows the prediction of the semiclassical diagonal 
approximation (23), which correctly describes the behaviour of the form factors for 
small r, but fails close to the Heisenberg time r = 1 as discussed at the end of Section 
[| in connection with the sum rule ([H]). The validity of this relation is demonstrated 
in Fig. |T0|b. It is interesting to note that all our data for systems with time-reversal 
invariance can be fitted very accurately (with A as a free parameter) using an ansatz 
which combines ([16|) and the semiclassical result (23) into 



K 0/1 (t) = K RMT (T)P 0/1 (Tt H ) 



(55) 
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Figure 11. Long-time behaviour of the group-element-specific form factors for the 
two-cell Sinai billiard, (a) shows i^ /i( T ) together with the prediction according to 
Eq. (^). The parameter (g) was obtained by a fit to the short-time behavior in the 
vicinity of the Heisenberg time r = 1 (see text), (b) shows the difference of the two 
form factors, which contains all essential information. 



We cannot further substantiate this expression analytically. 

Due to the poor outcome of the diagonal approximation in the present case, we 
cannot directly determine the mean resonance width from the decay constant by 
matching (23) and fl3"0|) at r = 1. Lacking a better semiclassical theory we fit (g) 
to our data, and we do so in the vicinity of the Heisenberg time. For the fitting we have 
actually used the fact that (|5^) as shown in Fig. [ll] represents our data up to and slightly 
beyond the Heisenberg time very well, (g) was determined as the value for which the 
long-time expression for the form factor matches smoothly to this ansatz. We prefer this 
procedure to a standard least-square fit over a large time interval, since it emphasizes 
that the value of (g) is at least implicitly contained in the short-time behaviour of the 
form factor. 

With (g) obtained in this way, fl3~0D describes K /i(t > 1) very well (Fig. [H]). Since 
the sum of the two form factors K + K\ is constant according to (|16D and Fig. [lO|b, all 
information is contained in the difference of the form factors which is shown in Fig. [lib. 
We regard the good agreement over a very long time as numerical evidence in favour of 
the presented theory for the long-time behavior of the form factor although it contains 
(g) as a fit parameter. 



6.3. Quantum graphs 

In this subsection we construct and investigate a two-cell system consisting of a 
quantized graph. It was recently shown [2TJ that quantum graphs exhibit the common 
quantum signatures of chaos and allow for a formally semiclassical description on the 
basis of non-deterministic classical dynamics. 

A graph is defined by v = 1, . . . , V vertices and 2B directed bonds connecting them. 
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The bond b with length is understood to lead from vertex v(b) to v(b), b being the 
reversed bond (Lg = Lb). On each bond we use a coordinate Xb with Xb = at v(b), 
Xb = Lb at v(b) and xi = Lb — Xb- The wave function 4> b (xb) satisfies the Schrodinger 
equation (H = 2m = 1) 

d -|2 \ 

— + k 2 ) <p b {x b ) = . (56) 
dx b \ J 

At the vertices, boundary conditions are chosen such that the current is conserved 
and the resulting Hamiltonian is self-adjoint and time-reversal invariant. Following the 



definitions in [21], we require (i), that the wave function is continuous across all vertices, 



i.e., it has the same value in all bonds b connected to some vertex v 

0fe(O) = 4>v=v(b) , (57) 
and (ii), that the sum of the momenta in these bonds vanishes 

E WJTT-MO) = 0. (58) 

b 



The eigenvalues of the so-defined graphs can easily be found numerically ]^TJ. For the 



unit cell, we have interconnected V = 10 vertices using B = 20 bonds such that each 
vertex is the intersection of exactly four bonds (inset of Fig. |T2|b). In this case the 
classical dynamics is particularly simple: On the bonds, there is free motion at speed 
2k and each vertex scatters the particle into any attached bond with equal probability 
1/4. On basis of such classical dynamics, a formally semiclassical quantization can be 
formulated which turns out to be exact in this model. 

One of the bonds connecting the two pentagonal layers in Fig. |T2"|b was sectioned, 



and both ends connected to a second identical unit cell, such that both cells form a 
ring with translation invariance in the direction "normal" to the pentagonal layers. The 
bond lengths Lb of the unit cell are chosen as random numbers, such that the reflection 
symmetry is broken. The total length is normalized according to J2b ^b=i = 17 so that 
the mean level spacing in k of the unit cell is unity. Therefore it is advantageous to 
use the wavenumber k and the path length L, instead of energy and time, as conjugate 
variables for the semiclassical description. The Heisenberg time is thus replaced by the 
Heisenberg length Lh = 2n and dimensionless time is introduced as r = L/Ln. In the 
ergodic regime, the escape rate from the unit cell (again with respect to unit path length 
instead of unit time) is simply given by the inverse total length of the graph, A = 1/tt. 

It is a particularly favourable feature of the model that neither this rate nor the 
number of quantum channels connecting the two cells depend on energy. Therefore, 
we can average over arbitrarily large energy intervals. We have computed the form 
factor from the 10,000 lowest doublets after dividing the spectrum into groups of 40 
doublets each. Fig. O compares the data as in the previous section to the diagonal 
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Figure 12. (a) Form factors K /i(t) for the two-cell quantum graph compared to the 
semiclassical diagonal approximation and the ansatz as in Fig. |To|a. (b) compares 
the difference of the form factors to the theory for the long-time behavior with the 
paramter (g) obtained as in Fig. [fl[ The insets show the analogue of Fig. p"l| a and the 
topology of the unit cell. 

approximation, the ansatz (|55|) and the long-time theory (f|l]). The results correspond 
to those for the Sinai billiard, but due the larger amount of data the agreement with 
our theoretical predictions is even closer. 

6.4- The quantum kicked rotor on a torus 

The kicked rotor belongs to the class of one-dimensional systems that are rendered 
classically chaotic only by a periodic driving. Its phase space is spanned by an angle 
and an angular-momentum variable and therefore has the topology of a cylinder. The 
nonlinearity of the potential is restricted to its time-dependent component and is 
controlled by a perturbation parameter. Accordingly, the classical dynamics crosses 
over smoothly from integrability to global chaos with increasing nonlinearity parameter, 
thereby following the KAM scenario |33|. The phase space of the kicked rotor 
is periodic also with respect to its non-cyclic coordinate, namely along the angular- 
momentum axis. 

Quantum-mechanically, the classical angular-momentum period coexists with h as 
a second independent action scale. If both are commensurable, then also the quantum 
kicked rotor is periodic with respect to angular momentum and can serve as a model for 
solid-state-like systems with discrete spatial translation invariance. Since in the periodic 
case, the cylindrical phase space may be regarded as being bent back to itself, this variant 
of the model is referred to as the 'kicked rotor on a torus'. It is this case which we shall 
discuss below. If the two angular-momentum periods are incommensurate, the quantum 
eigenstates are generally localized. In this case, the kicked rotor provides a model for 
Anderson localization in disordered systems 0, We will not consider it here. 
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The kicked rotor is defined by its Hamiltonian 

H(l,-&-t) = ^-—L+V a)k (-&) £ <J(t-mr). (59) 

m=— 00 

As a consequence of the periodic time dependence, spectrum and eigenstates are 
adequately discussed in terms of quasienergies and Floquet states, respectively In 
addition, the kicked rotor may possess two independent twofold antiunitary symmetries, 
both resembling time-reversal invariance. In order to break them in a controlled manner, 



an angular-momentum shift A has been introduced, and the potential is chosen as j34 

= k 



IT 1 7T 

cos(a— ) cos & + - sin(a— ) sin 2§ 



(60) 



Here, the global prefactor k determines the degree of nonlinearity. The appropriate 
classical measure of nonlinearity, however, is the parameter K = kr. If K ^> 1, chaotic 
motion prevails and angular momentum diffuses without restriction by KAM tori. 

The ratio of the classical to the quantum period of action is determined by the 
parameter t/Att. If it is rational, i.e., if r = 4np/q, with p, q coprime, a unit cell along 
the angular-momentum axis arises that accomodates q quanta of angular momentum. 
The number of quasienergy levels per unit cell is then also q. We set p = 1 and, in order 
to avoid an unwanted symmetry of the unit cell, require q to be odd. 

According to Bloch's theorem, the spatial periodicity implies the existence of an 
additional constant of the motion, the Bloch phase 9. It appears explicitly in the 
symmetry-projected Floquet operator P2"] . 

(Z' |*7(0)|O = exp^-27ri^(/-A) 2 ^ 

1 9-1 

x _ e -iV a , k ([e+27Tn]/q) e i(l-l')(e+27Tn)/q_ ^j-Q 
1 n=0 

A restriction of the lattice to a finite number of A" unit cells, with cyclic boundary 
conditions at the ends, amounts to discretizing the Brillouin zone so that it comprises A" 
equidistant values 6 m = 2irm/N, m = 0, . . ., N—l, of the Bloch phase. The independent 
parameter A^ corresponds to the number of levels per band, the total number of levels 
in the spectrum is Nq. 

A system with two unit cells is constructed simply by setting N = 2. In contrast to 
the systems discussed above, the resulting model does not possess a bottleneck between 
its two compartments, neither in configuration space nor in phase space. A reduced 
exchange between them therefore comes about solely by slow diffusion. The exchange 
rate A is determined by the diffusion constant D = k 2 /2 (valid if K ^> 1) through the 
simple relation A = D / (2a 2 ) derived in Appendix A . In the quantum kicked rotor on a 
torus, the integer q represents the dimensionless size of the unit-cell, to be substituted 
for a. 
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The most interesting parameter regime to be studied numerically would be one where 
Kq{t) exhibits a positive peak around r = 1, the feature indicating quasidegeneracy in 
the spectrum. For this peak to emerge, the exchange between the cells must be slow. 
Since there is no bottleneck in the kicked rotor, this can only be achieved through a 
small diffusion constant. More precisely, it requires that the Heisenberg time should be 
small against the Thouless time. Measured in units of the discrete time steps of the 
kicked rotor they are, respectively, % = q and n D = N 2 q 2 / (irD), so that the condition 
for quasidegeneracy to occur reads 

N » W (62) 

At the same time, it should be avoided that localization becomes effective even within 
the unit cells, in order to separate the signature of classical diffusion from the direct 
quantal effect of disorder in the spectrum. The localization length should therefore be 
kept large compared to the size of the unit cell, 



f « — > q. 



(63) 



Clearly, both conditions, (p2[) and (|63|) , can hardly be met simultaneously if N is fixed 
and small. With N = 2, little freedom remains since, in addition, being close to the 
classical limit and well within the classically chaotic regime requires both q and k to 
be large. We found that q = 45 and k — 10 represents an acceptable compromize. 
The resulting diffusion constant, corrected for oscillations occurring if K^l [32], is 
D = 23.23. We substitute the Thouless time td = n^/q = 2.467 for the time constant 
1/(A£h) of the exponential equilibration. 




Figure 13. Group-element-specific form factors (a), as described by Eqs. (23) and 
(|52"|), and corresponding cluster functions (b), for the quantum kicked rotor on a torus, 
compared with the theory according to Eqs. (23) and ( |52"| ) (dashed). The parameter 
values are q — 45 and k = 10, corresponding to an exchange rate A = 0.258. 



In Fig. [13|, we compare the form factors Kq/i(t) (panel a) and corresponding cluster 
functions y /i( r ) (b) obtained for the quantum kicked rotor on a torus with parameters 
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as above, with our theory, Eqs. (23) and (|52|). For the evaluation of the theory, we have 
used the relations cited to determine the decay rate directly from q and k. No fitting 
was involved. The data cover ten Heisenberg times and thus reach far into the quantum 
long-time regime. In the form factors, we can clearly discern the three time domains 
discussed, the initial phase of chaotic diffusion where Kq(t) is strongly enhanced while 
Ki(t) remains close to zero, the sharp positive peak of Kq(t), reaching almost twice 
the asymptotic value, and the saturation regime where K (t) and Ki(t) approach their 
common asymptote from above and below, respectively. The cluster function for r ^ 1 
represents the regime of long times or small splittings in a different manner. Both 
plots give evidence that the theory provides a quantitative description of the two-point 
correlations over all time/energy scales. 



6.5. A random-matrix model 

In the present subsection, we consider a simple model for a reflection-symmetric 
double well system with chaotic dynamics, with and without time-reversal symmetry, 
constructed in the spirit of Refs. [|11], |I2"[| . It will be shown that in certain cases, this 
model qualitatively reproduces the features of the form factor discussed in Sections f| 
and [|. 

We consider a Hamiltonian of the form 



H= ° I (64) 




where H represents the internal dynamics of either cell in an iV-dimensional Hilbert 
space, and V their coupling via M N channels of the connecting duct. Note that in 



contrast to the models discussed in || [LI], [12] , we require the two blocks on the diagonal 
to be identical. 

We model Hq as an Nx iV random matrix distributed according to Dyson's Gaussian 
ensembles, P(H ) dH oc exp(— Tr Hq/A) dH . It is assumed that iV — > oo. The N x N 
matrix V has the form 

N vA 

V k i = 5 kl — — forfc = l,...,M (65) 

M 7l z 

and zero for k > M. Here, M -C N is the number of matrix elements coupling the two 



wells, v parameterizes their strength, and A is the mean level spacing of H , A = -n^J (3/N 
with (5 = 1 in the Gaussian Orthogonal Ensemble (GOE) and (3 = 2 in the Gaussian 
Unitary Ensemble (GUE). 

The Hamiltonian H has a twofold symmetry. Its eigenvalues can be classified 
according to parity p and appear as doublets r a ^ u with a = 1, . . . , N and v — ±. 
According to Eq. we write r a> ± = R a ± r a . The form factor is then given by (cf. 
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Eq. m 



1 / 1 N \ 
Ko/l (r) = - ± ( E cos(47rrr a ) \ . (66) 

For large times (r ^> 1), ifo/iCr) m &y be calculated by evaluating the doublet splitting 
2r a within degenerate perturbation theory. Denoting the eigenfunctions of H by <ft a 
(with components (p au ), one has (with (did) = A -1 ) 

(67) 

Substituting (|67[) into (|66"D, it remains to average over the eigenfunctions (fi a of H . The 
statistical properties of the eigenfunctions <p a depend on the ensemble considered. In 
the GUE, the amplitude u = N\<fi au \ 2 is distributed according to P(u) = exp(-u). In 
the GOE, the corresponding distribution is P(u) = (2ttu)~ 1 ^ 2 exp(— u/2). 

We first consider the case M = 1, where the two wells are coupled via a single 
matrix element. For r ^> 1, one obtains for Kq/i(t) 

1 1 f 2 (vt/tt)- 1/2 for /3 = 1 

K 0/1 r )»i±-J 4V ^ (68) 

2 2 I 1/[1 + (4-yr/Tr) 2 ] for/5 = 2. 

For /? = 1, this expression reproduces the long-time r~ 1//2 decay of Eq. fl44[) . For (3 = 2, 
Eq. ( |68D reproduces the r~ 2 decay for large r implied by Eq. ([47D. 

For p = 2, the matching procedure discussed in Section |^ yields an analytical 
expression for -Ko( r ) valid for all time scales. In Fig. |14|a, we compare this expression 
(Eqs. (23), (|52|)) with results of simulations of the model (0). Shown is K (r) as a 
function of r for an ensemble of random matrices with N = 80, v = 0.5, 1 and 2 in the 
GUE (full lines), as well as Eqs. (23) and (0). The constant v may be determined from 
(g) by comparison of Eq. flBlf) with Eq. fl52|), where c = 8ir(g). 

We find good agreement between the results of the simulations and Eqs. (23), ([52]) . 
We have, however, not attempted to evaluate the small-r behaviour of the form factor 
for the model QS^D analytically. For large r, on the other hand, it is clear that Eq. 
is a good model for the form factor: as pointed out above, Eqs. flBS|) and (|5*^ ) coincide 
for large r. 

In the present model, it is also possible to consider larger M, 1 < M N. In this 
case, for (3 = 2, the quantity 

N M 

M = T?EI^| 2 ( 69 ) 

m v=\ 

is distributed according to 
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Figure 14. (a) Form factor Kq(t) in the GUE for three different values of the coupling 
strength (v = 0.5, 1, and 2), and for N = 80 and M = 1 (full lines). Also shown are 
the theoretical results, Eqs. (23) and ( |52| ) (dashed), (b) Form factor Kq(t) — 1/2 in 
the GUE for v = 0.2, TV = 80 and M = 1, 2, 6, and 10 (full lines), compared to the 
asymptotic theoretical result, Eq. (71) (dashed). 

We require M <C N since for M = N, one has u = 1 due to normalization of the wave 
functions. In the GUE, the form factor is then given by 

1 1 r°° /4v \ 
K (t)^- + - du cos — ru )P(u, M) 

2 2 JO V 7T ' 

1 1 / S4VT\2\ ~ M I 2 ( /AVT\\ 

= - + - 1 + ( ) cos Matan . 71 

For M — 1, Eq. ((68|) is reproduced. For 1 < M < Af, the form factor does not decay 
monotonically for r > 1 but exhibits oscillatory behaviour (Fig. |T4|b). 

7. Conclusion 

Chaotic systems with two weakly connected cells, elementary as this concept may 
appear, form a paradigm for a large class of physical situations and exhibit a surprisingly 
rich behaviour. In this paper, we have shown that it is determined essentially by two 
parameters. One of them can be identified with the time required for the respective 
populations of the cells to equilibrate. It specifies the position between the extreme of 
a large opening that hardly restricts the exchange, and the opposite one of two almost 
isolated single cells. The second relevant parameter is a measure of the difference in 
shape between the cells, ranging from exact symmetry to its complete absence. 

On basis of the results obtained in this paper, we can draw a clear picture of the 
spectral two-point correlations in this two-dimensional parameter space. In the case of 
an effective communication between the cells, the presence or absence of symmetry is of 
little relevance for the spectral statistics. It is then only the slight retardation of ergodic 
coverage that becomes manifest in the level correlations. The result is a reduction in the 
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area enclosed by the initial minimum (correlation hole) of the form factor, indicating 
an increase of randomness in the spectrum which can be completely accounted for by 
semiclassical considerations ||. 

The case of two almost isolated cells lacking all symmetry can be trivially understood 
from a random-matrix point of view. We are then dealing with the superposition of 
two spectra that are nearly mutually independent but exhibit the same statistics. Here, 
random- matrix theory simply predicts a doubling of the time argument of the form factor 
P5| ] , in agreement with the semiclassical approximation in the limit of slow equilibration. 

If, in contrast, the two cells are symmetric, the formation of doublets introduces an 
additional feature in the spectrum. The corresponding positive correlations are reflected 
in the form factor as a maximum in the vicinity of the Heisenberg time. In the limit 
of long exchange time, the form factor at this maximum reaches twice the asymptotic 
value to which it decays subsequently from above, relative to its value at t — 0. In the 
case of exact symmetry, we can quantitatively account for this peak in the standard 
form factor. Simultaneously, there is a depression in an analogous statistic that refers 
to transport from one cell to the other, rather than to return to the initial one. 

The crossover from full to completely broken symmetry, as a function of some 
symmetry-breaking parameter, can be included in the semiclassical theory if a few 
plausible additional assumptions are made. In accordance with corresponding work 
on spatially periodic systems with slight disorder ||3~5] , this approach implies that the 
peak in the form factor should decay exponentially both with the typical difference in 
action between symmetry-related periodic orbits in the respective cells, and with time. 
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Appendix A. Discrete diffusion 

The picture of a two-cell system with the topology of a ring suggests to unroll the ring 
so that an infinite chain is formed. If the two cells are translation symmetric, then each 
of them represents a unit cell of this periodic lattice, otherwise the unit cell comprises 
both cells of the ring. In this extended topology, the condition of slow exchange between 
the cells implies that the picture of homogeneous diffusion breaks down on the scale of 
the lattice constant. We are therefore in a regime opposite to that considered in Ref. KJ. 
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The spreading along the chain is now determined by a master equation for the 
probability to be at site n, instead of a diffusion equation for the probability density 



P n (t) = XP n -!(t) - 2XP n (t) + AP n+1 (t), 
In a more concise notation, it reads 

(A-AA 2 )p(t H 0. 

Here, P(t) denotes the entire infinite vector of the P n (t), and 



(Al) 
(A2) 



\ 



. . 1 
1 -2 1 

1 -2 1 

1 -2 1 
1 . 



V 



(A3) 



is the discrete Laplace operator. The lattice plane waves 

1 



4>r, 



^2-Kinm/N 



N 



m = 0, ...,N- 1, 



(A4) 



where for the sake of normalizability we have reintroduced cyclic boundary conditions 
with a period of N chain elements, solve the stationary eigenvalue equation 

with eigenvalues 

2nm 



7m = 2 (cos v 
For a localized initial state 

P n {0) = S nmo dN 



(A5) 
(A6) 



-i N-l 



(AT) 



m=0 



the time evolution reads 



Pn(t) = 



1 



N-l 



„ E ^,n« 
7V m=0 



,7m* 



JV-1 



m=0 



mn 



-Eex P (2 7 ri— + 2A 



AT 



27rm 

COS— -1 



(A8) 



By Poisson resummation, this becomes 



p n(*) = -5JF E / di/exp 27ri 

oo 

3_ * E I|n+mJV|(2A£), 



m + 



iV 



i/ + 2A 



2vrz/ 

COS— -1 



(A9) 
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where l n (z) denotes the modified Bessel function of integer order n JHJ. The spreading 
over the lattice, as described by P n (i), represents a discrete diffusion process. If we go 
to the continuum limit by defining x = na, L = Na, and letting the lattice constant 
a — > 0, we recover continuous diffusion with periodic boundary conditions, 

/ s Um 1 ^ ( (x + mL) 2 \ /A , 

p(x,t =-P B t £ «P - on. ■ Al ° 

The diffusion constant is D = 2Xa 2 . In performing the limit, we have used the 
asymptotic form of the l n (z) for large argument z |37] and expanded it for large order 



n. The two-cell solution, Eq. (0), is retained by setting N = 2 in Eq. (A8). 
Appendix B. Doublets as discretized bands 

As on the level of the classical dynamics, it is instructive to consider also the quantum 
two-cell system as the unit cell of an infinite chain. From this point of view, the doublets 
r Q) ± come about by discretizing continuous bands to a "Brillouin zone" with only two 
points. The simplest possible interpolation between these points assumes cosine-shaped 
bands, 

r a(/-0 — Ra + r a COs(7T/i), \l = 0, 1. (Bll) 



Equation (|BII| ) can be justified by the fact that it imposes no more information on 
the shape of the bands than is available, namely their first two Fourier coefficients. 
Cosine-shaped bands result also from diagonalizing a tight-binding Hamiltonian with 
translation invariance. For two sites this is 

Rn,n' = Ra5( n -ri) mod2 

+ 2 r "(^(™-«'-l)mod2 + ^(n-n'+l)mod2)) 

n,n' = 0,l. (B12) 

We have defined the parameters of this Hamiltonian in such a way that Eq. (pil|) gives 
its eigenenergies. Inserting them in Eq. (27) and performing a Poisson resummation 
results in 

oo -/V d 

«n(r)= E E e_27riT ^ in " 2mj 2--n(2vrrr a ), (BI3) 

m=— oo a=l 

where Jfc(z) denotes the ordinary Bessel function of order k. We introduce a diagonal 
approximation with respect to the band index a, as in Section |5], and obtain the 
corresponding form factors as 



I 



Kn{r) = — E 

7V d a =l 



E (-l) m h m -n(2nrr c 



m=—oo 



2 



(B14) 
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Invoking the sum rules J2kL-<x>(~ l) fc J 2 fc(-2) = cos z, Y,T=~oo(~^) k ^2k-i(z) = sinz |38[ 
we recover Eq. fl29|). 



Appendix C. Quantization of the Z-shaped billiard 

Quantization of a billiard amounts to solving the Helmholtz equation 
/ d 2 d 2 



\o,* + W + k )^ x > y) = ° (C15) 

with Dirichlet boundary conditions on the billiard circumference and a dispersion 
k 2 = 2mE/h 2 . In this appendix we describe a specific quantization method for the Z- 



shaped billiard discussed in Section |0]. It is based on the scattering approach [g6j. [2(J • 
Consider a subdivision of the closed double billiard into two open halves (Fig. 0). Each 
of them represents a chaotic scatterer attached to the end of a semi-infinite waveguide. 
Within the waveguide, quantization of transverse momentum, k y>n = nir, n = 0, ±1, 
±2, . . ., implies that there are N = [k/n] ([. . .] denoting integer part) open channels 
with real longitudinal momentum k x<n = {k 2 — k 2 ^) 1 / 2 , such that the two scatterers are 
described by iV x iV scattering matrices S l and S T , respectively. 

The secular equation for the eigenvalues of the full billiard then reads 

det(I - S\k)S r {k)) = 0. (C16) 

In order to construct the S l ^ T for the billiard halves 1 12611, we start from the 2N x 2N 



transfer matrix for a quarter Sinai billiard open on both sides @,|], 

(C17) 



rt 1 r — t rt 1 

rV t- 1 



Here, t and r denote the N x N matrices of transmission and reflection amplitudes, 
respectively. Due to the spatial reflection symmetry with respect to the diagonal, the 
two entrances of the billiard are equivalent. 

The transfer matrix T w for a waveguide of length L/2 consists of phase factors 
exp (±ik x<n L/2) along the diagonal. The letter- Z-like fashion in which the two halves are 
assembled is accounted for by a third factor T z with appropriate phases ±1 0, [| along 
its diagonal. It is included in the transfer matrix for one of the sides, e.g., T 1 = T e T w , 

The scattering matrices for the billiard halves closed on one side are obtained from 
T l and T r by requiring incoming and outgoing amplitudes to cancel across the openings 
where Dirichlet boundary conditions are to be enforced, 

(C18) 
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Here, A and —A refer to the amplitudes at the ends to be closed, and B and C to the 
amplitudes on the opposite sides. The latter are related by C = S l ^ T B, invoking the S 
matrices sought for. Solving for them, one finds 

S l/ * = — (J12 + T l £)-\T^ + T l 2 { r ). (C19) 

l/r 

In obvious notation, , i,j = 1, 2, refer to the four NxN subblocks of the respective 
transfer matrices. 




Figure Al. The cumulated eigenvalue density of the double billiard (steps) compared 
to the Brownell formula (dashed line). The geometric parameter values are R = 10 
and L = 10. 



In order to check the quality of the quantization procedure, we compare, in Fig. [Ad 
the numerical result for the cumulated eigenvalue density with the Brownell formula |29[ 



for wavenumbers in the interval tc < k < 2tt, and R = 10. The agreement is satisfactory 
up to k m 5. For larger wavenumbers, quasidegenerate pairs of zeros occur in the secular 
function with too small spacing to be resolved by the numerical procedure. We have 
therefore discarded data with k > 5. To achieve better statistics in the evaluation of 
spectral correlations, we have varied L within an interval AL amounting to a few percent 
of L. 
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